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Motivated by recent studies of the one-bubble inflationary universe scenario that predicts a low 
density, negative curvature universe, we investigate the Euclidean vacuum mode functions of a scalar 
field in a spatially open chart of de Sitter space which is foliated by hyperbolic time slices. When we 
' consider the possibility of an open inflationary universe, we are faced with the problem of the initial 

0^ condition for the quantum fluctuations of the inflaton fleld, because the inflationary era should not 

last too long to lose every information of the initial condition. In the one-bubble scenario in which an 
open universe is created in an exponentially expanding false vacuum universe triggered by quantum 
O ■ decay of false vacuum, it seems natural that the initial state is the de Sitter-invariant Euclidean 

vacuum. Here we present explicit expressions for the Euclidean vacuum mode functions in the open 
chart for a scalar field with arbitrary mass and curvature coupling. 
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I. INTRODUCTION 
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Recently there appeared not a few observations which suggest our universe has negative curvature, i.e., fig ~ 0.1 
101 . Accordingly, papers comparing theoretical predictions of open universe models with observational data such 
as COBE Q have been appearing in the context of inflationary universe models in cosmological models with 
topological defects or in a general context by assuming a power-law type primordial density perturbation spectrum 
1^. However, in the standard inflationary universe paradigm, it is generally believed that an open FRW universe with 
small perturbations is hard to be realized in a consistent manner H . 
h"*' One possible consistent scenario is the creation of an open universe from an exponentially expanding false vacuum 
rr| dominated universe |0-|lO|. Originally this idea was proposed by Gott |l|]. In the standard inflationary universe 
^ • scenario, the horizon problem is solved by a large amount of expansion of space. A homogeneous patch initially of a 
horizon size expands exponentially and our present horizon size will be inside such a homogeneous patch. However 
in this context, the flatness (or entropy) problem is solved at the same time when the horizon problem is solved. 
■ Therefore the spatial curvature at present time is inevitably decreased by the expansion of universe. This is the 
problem of the standard scenario of inflation if we attempt to construct an open universe model with 1 — f^o = 0(1). 
On the other hand, if we consider a bubble nucleation in the sea of false vacuum which is described by de Sitter space, 
the interior of a bubble has the 0(3, 1) invariance owing to the 0(4) symmetry of the Euclidean bounce solution which 
represents the tunneling process |^ . Thus the horizon problem is automatically solved. Though the spacetime has 
the exact 0(3, 1) invariance in the lowest order description, quantum fluctuations around the classical background 
may give rise to cosmological density perturbations to explain the large scale structure of the universe. However if 
the vacuum energy never becomes a dominant component of the cosmic energy density after nucleation, the universe 
will be curvature dominated from the beginning and never recover to a hot FRW universe. Thus the entropy problem 
cannot be solved. To solve this problem, a secondary inflation in the bubble is required. The essential difference of 
this scenario from the standard scenario is that a horizon size patch whose size is approximately equal to the curvature 
scale at the onset of the secondary inflation may not become much larger than the present horizon scale. In such a 
case, we will have a sufficiently homogeneous open universe at present. This implies that memories of the quantum 
state of the universe at the beginning of the secondary inflation will not be erased but will directly affect the observed 
large scale temperature and density fluctuations p-[lO||. 
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Therefore the quantum state of a field, 4>, inside the nucleated bubble, especially that of the inflaton field of the 
secondary inflation, is to be examined. A pioneering study of this subject was done by Rubakov ||l3|] and a formalism 
which respects the 0(4) symmetry of the tunneling background was developed by Vachaspati and Vilenkin [Q. 
Meanwhile we developed a formalism based on the multi-dimensional tunneling wave function ||l5[] and applied it to 
the 0(4)-symmetric bubble nucleation without gravitational effects and with gravitational effects However, 
all of these previous works remained in a rather formal level in the sense that techniques to investigate the quantum 
state which have practical applicability to a general situation have not been developed. Recently, we have succeeded 
in giving a general and practical method to obtain the quantum state inside a nucleated bubble in flat space fl^ and 
its application to the case of the quantum state of the tunneling field itself is in progress ||l9| . Thus an extension of 
this method to the case when gravity comes into play is now to be formulated. 

As a first step, we consider the case when the gravitational back reaction effect can be neglected so that the 
background metric may be fixed to that of de Sitter space. Then the spacetime inside a bubble is described by a 
spatially open chart of the de Sitter space and the inflaton field will take a constant value on a hypersurface of the 
spatially open time slicing. Further, if the tunneling field which causes the bubble nucleation has no interaction with 
(/), the quantum state of (j) will not be affected by the tunneling process at all. Then, provided that the quantum 
state of before tunneling is in the Euclidean vacuum, which should be a good approximation if the preceding false 
vacuum inflation lasted long enough, the quantum state inside a nucleated bubble will remain so. Hence, it is required 
to describe the quantum state by the mode functions in an open chart of de Sitter space. Once we know a method to 
obtain the Euclidean vacuum mode functions in the open chart, it should be fairly straightforward to extend it to the 
case in which the mass of (j) changes in time due to its coupling to the tunneling field or when the geometry deviates 
from the exact de Sitter space at later stages. 

In a spatially flat or closed chart, the field operator can be relatively easily decomposed into the spatial harmonics 
and the mode functions corresponding to the Euclidean vacuum are well known [ pO[ . However, the Euclidean vacuum 
mode functions in a spatially open chart have not been known except for the massless conformally coupled case 
Here we give an explicit expression of vacuum mode functions for a scalar field with arbitrary mass and curvature 
coupling, including the massless minimal coupling limit. 

The paper is organized as follows. In Section II, we quantize a scalar field on the hyperbolic time slices which 
foliate two distinct open charts of de Sitter space (see Fig. 1). We define a vacuum state there by requiring that the 
positive frequency functions be regular on a hemisphere of the Euclidean de Sitter space where the complexified time 
coordinate is negative pure imaginary. Then we derive an expression for the Wightman function for this vacuum state 
in the series form. In Section III, by analyzing the behavior of the Wightman function for thus obtained vacuum 
state, we show that it is in fact the de Sitter-invariant Euclidean vacuum, the so-called Bunch-Davies vacuum [po[ , 
provided that the effective mass of the scalar field is greater than a critical value which corresponds to the conformally 
coupled massless case. However, we also find that our expression for the Wightman function does not coincide with the 
Euclidean vacuum one for mass smaller than the critical value. In Section IV, by carefully analyzing this discrepancy 
in the small mass case, we find that there exist a set of modes which have finite Klein-Gordon norms on regular 
Cauchy surfaces, say on spatially closed time slices, but which cannot be quantized on the open charts because of 
the divergent Klein-Gordon norms on the hyperbolic time slices. Then normalizing these modes on a spatially closed 
hypersurface and analytically continuing them back to the open charts, we obtain a complete description of a scalar 
field and the Euclidean vacuum there in terms of the orthonormalized positive frequency functions, irrespective of 
its mass. In Section V, we consider the special cases of a massless conformal scalar and a massless minimal scalar 
in which the series expression for the Wightman function can be summed up with elementary algebra. The results 
are found to be in perfect agreement with our analysis for general mass in the preceding sections. In Section VI, 
we summarize our results and discuss their implications. Finally, Appendix A explicitly evaluates the Klein-Gordon 
norms of the mode functions on a closed slice for which their norms on a hyperbolic slice are finite, to show their 
equivalence and Appendix B gives a proof of a mathematical formula which plays a central role in the analysis of the 
set of modes whose Klein- Gordon norms diverge on a hyperbolic slice. 

II. SCALAR FIELD ON OPEN DE SITTER SPACE 

To begin with, we introduce a coordinate system which covers the whole Euclidean de Sitter space. The four- 
dimensional Euclidean de Sitter space is a four-sphere and can be embedded in the five-dimensional Euclidean space, 
(x^jX*), as a hypersurface which satisfies = {x^Y + Tlii^^Y ^ where is the Hubble radius of the de Sitter 
space. A fundamental coordinate system we use in the Euclidean region is defined as 
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Then the metric is represented as 
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Coordinate systems in the Lorentzian region are obtained by the analytic continuation, a;' 



Performing this, 



we find that the Lorentzian region may be divided into three parts, which we call R, C, and L, whose coordinates are 
related to the hmdamental ones by the relations. 
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and their metrics are given respectively by 



ds^c 



dsl = H-^ 



sinh^ tLi{dr]i + sinh^ rjidfl^)^ , 
+ cos^ tc{~drl, + cosh^ rcdft'^)) , 
-dtl + sinh^ tL{drl + sinh^ rLd9?)) . 



(2.4) 
(2.5) 



In the above, the regions R and L are connected through r € (— 7r/2,7r/2) along a path over the hemisphere 5a;° < 
(x*' > 0), say along p = 0, and the region C is joined to this hemisphere at p = 7r/2. Note that one may connect 
the regions R and L in a different way by passing through the other hemisphere '^x^ > (5° < 0), say along p — it, 
but we take the former choice given by Eq. ( |2.3| ) for later convenience. A schematic picture of how these coordinates 
cover the spacetime is shown in a conformal diagram in Fig. 1. The regions R and L covered by the coordinates 
{tR, tb) and (<l, ^l), respectively, are the two distinct spatially open charts of de Sitter space. We see that the whole 
Euclidean and Lorentzian de Sitter space is covered by complexifying the variables r and p. 

Let us consider a free scalar field in de Sitter space with curvature coupling constant ^. We expand the field operator 

as 



) = ^ \ atvUA{x) + d\uA{x] 



(2.6) 



where a bar denotes the complex conjugate and {uA(a;)} forms a complete set of mode functions labeled by certain 
indices A which satisfy the field equation, 



[g'^'V^V, - Mlfj] ua{x) = 0; M^fj := + 12(H\ 



(2.7) 



and normalized with respect to the Klein-Gordon inner product. The condition, oaIO) — for any A, determines the 
vacuum state associated with a specified set of mode functions. 

To find the Euclidean vacuum mode functions on the open chart, we write down the field equation ( ^.7| ) in terms 
of the coordinates in either of the regions R or L. Since these two regions are completely symmetric we have 
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where (t, r) = {tji^ra) or {tL^rL), 
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UA(i,r,r!) = 0, 



(2.8) 



a{t) = H~^sinhi, 
I d 1 



sinh^ r dr 
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and Lq is the usual Laplacian operator on the unit two-sphere. We write the eigenvalue equation for the operator 
— on a unit three-dimensional hyperboloid in the form [p2| , 

- L%irn{r,n) = {l+p^)Yp,rair,n). (2.10) 

The eigenfunction Ypim which is regular at r — is given by 



n) = fpi{r)Yi^{n) 

r{ip + i + i) p p-i-i/2 

T{ip+1) Vsinh 



Mr) - ^^^^:^^^C7/(cosh.) 



= i-iyj- sinh' r / , f ^) . (2.11) 

^ ' \l t: T{-ip + l + l) d{coshry \smhr J ^ ' 

where Yim{^) is the normalized spherical harmonic function on the unit two-sphere, T{z) is the Gamma function and 
Pj^iz) is the associated Legendre function of the first kind The eigenf unctions Ypim with real positive values of p 
form an orthonormal complete set for square-intcgrable functions on the unit three-hyperboloid. They are normalized 
as 







drsinh r / dV,Ypijn{r,^^)Yp,i,jn'ir,fl)^S{p~p')SwSmni'- (2-12) 



We now take the harmonic expansion of the mode functions, 

Uplm{t,r,il) = -^Xplm{t)Ypijn{r,il), (2-13) 

a[t) 

and look for a complete set of positive frequency functions {Xpim(t)} for the Euclidean vacuum. 

As well-known, the positive frequency functions of the Euclidean vacuum are characterized by their regularity on the 
3a;° < hemisphere of the Euclidean de Sitter space. To find such functions, let us first consider the mode functions 
in the region R. The general solution is expressed as a linear combination of the associated Legendre functions 
P^J'_i/2i^R) ^^'^ -^J^i/2('^^j)' '^hsre ze = cosht^. To find out their coefficients, we first ignore the normalization and 
consider a solution, 

) = P'J:{zr)- (2.14) 

where we have introduced v' := v—1/2 for notational simplicity. Note that this function would be a natural candidate 
for the positive frequency function if the region R were the whole universe. Imposing the above regularity condition 
on Xp^\ we find that the analytic continuation of it to the region L amounts to going through the branch cut [—1,1] 
of the associated Legendre function P^{z) from 3z < to 3z > and then to z < — 1 on the real axis. Thus in the 
region L we have 

^(R) ^ e--^'Pj(-Zi) 

^ ^.(-,W)p^P(,^) _ ^£!!il(!^±i^Q^P(,^), (2.15) 

TT 

where = cosht^ and Qf-{z) is the associated Legendre function of the second kind [^. Replacing the role of R 
and L, we obtain Xp"^ in same way. With this choice of the mode functions, their regularity on the 'ix^ < 
hemisphere is automatically guaranteed because the harmonic function Ypim is regular at r^i^ = p — 0, hence 
regular for < p < ■n/2. 

As the mode functions we have obtained are not yet normalized, the Klein-Gordon inner products of these mode 
functions should be calculated to normalize them. The evaluation of the Klein-Gordon inner products must be 
performed on a Cauchy surface. As the choice of a surface is arbitrary as long as it is a Cauchy surface, we choose it 
in the following way. It consists of the three parts, (I), (II) and (III), where (I) is the rj^ < Vmax part oitfj — constant 
hypersurface for a large rmax, (II) is the one in which R is replaced by L, and (III) is a bridge connecting these two 
isolated parts. A schematic picture of this Cauchy surface is drawn in Fig. 2. If the contribution from the integral 
over the surface (III) can be neglected for sufficiently large Vmax, the Klein-Gordon inner product reduces to the 
summation of integrals over the surfaces (I) and (II). In such a case we have 
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(2a|)v„„(,,o),M|)yi„„..(r,o)) 



[{,4 



S{p - p')Sii'S n 



(2.16) 



As will be discussed in Section IV, this is not always correct but we assume so for the time being, 
in inner products of Xp^'' ^-nd Xp^^ 



Then the Klein-Gordon inner products of Xp^'' ^-nd Xp^^ ^I'e calculated to be 



= ((X^'^^X^'^^)) = ^e-^'(cosh27rp-cos27r;.'), 

where we have used the fact that 



(2.17) 



(2.18) 



and the Wronskian relations among P^^^ and Q^,^^ |2^. In the above and in the rest of this section, we assume 
to be real, i.e., —1/2 < i/' < 1 (0 < v < 3/2), in order to avoid inessential complexity. Extension to the case z/' 
is imaginary, i.e., ^v' = —1/2 {v = pure imaginary) is straightforward. It then turns out that the following linear 

combinations of Xp^"* ^-nd Xp"^ ^Iso mutually orthogonal: 



Xp,+ 



Xp + Xp 



(1 + e-^P+'''''^)T{v' + l + ip) 



1 



T{v' + l + ip) 
1 

[ T{v' + l + ip) 



^JM + -(i-e^"-^^")g:.'^(zfl) 

TT 



^J(^L) + -(l-ef— -)Q';(zi) 

TT 



Xp,- 



(«) (i) 
Xp - Xp 



(1 - e-^P+'''-^'-)T{v' + 1 + ip) 



1 



If we note the relation. 



T{v' + l + ip) 
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[ ^V{v' + l + ip) 



TT 



TT 



(2.19) 



2i sinh TTp 



pip _ ^{t^' + l + ip) p-tp 
T{iy' + l-ip) 



(2.20) 



we may re-express Xp.o- in terms of P^f alone as 



Xp,. 



o-n'P 



-7rp 



2 sinh TTp \T{iy' + ip+l] 



r(i/' -ip + i) 



(2.21) 



2sinh7rp \T{v' + ip+l) 



7-,P'Ji^L) - % , '"\., P;''{zl) 
1) 1 U^' — ijD -I- 1 I 



Note that these mode functions have the symmetry X-p.o- = Xp.o- (c = i), hence are even functions of p. Their 
Klein-Gordon norms are evaluated to give 



4 (coshTrp — crcosi^'Tr) 

TT \T{l'' + l+ ip)\ 



(2.22) 
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From these results, the field operator is now expanded as 

(f){x) ^ J dp ^ (apahnVpahn{x) + a^p^i^Vpalmix)^ , 



(2.23) 



{T,l.7n 



where and in what follows, we use the symbol v\ to denote the orthonormalized mode functions. In the present case, 
they are given by 



(2.24) 



with (i, r) being either (t/j, rn) or (t^, r^). In what follows, we suppress the subscript R (or L) for notational simplicity 
unless ambiguity arises. 

Then, for both x' and x" in the region R (or L), the Wightman function is given by 

poo 1 

a{t')a{nG+{x\x")^ / dp ^ — XpAt')ypim{r' ,n')xpAt")yptm{r" , 



where we have used the completeness relation for Ypim, 



(2.25) 
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Lra 



p sinpC 
27r2 sinh C 



cosh^' := cosh r' cosh r" — sinh r' sinh r" cos 9, 



(2.26) 



with COS0 being the directional cosine between fl' and Q" . Using the relations (2.18) and (2.2C), and the symmetry 
X-p,cr — Xp.cr again, the above integral is written explicitly as 
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sva{ip + v' )t: 
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sm(ip — V )t: 



2 sini^'TT 



TT sin(zp + v')!: 
2 sin v'-K 



TT sm(ip — j/')7r 



— e'^^'p;;^^(z')g;'^(z") 



(2.27) 



In deri ving the above expression, we have assumed j/ = i^' + 1/2 to be real. However, it can be shown that the final 
result (2.27) equally holds for imaginary v without change. Note that this is true irrespective of the sign of v, i.e., 
for v — zLi\iy\. Note also that p = is no longer a pole for i'' T^integer, which we assume in the rest of this section. 
The cases i/' ^ {i/ ~ 1/2; massless conformal) and ly' = 1 {v = 3/2; massless minimal) will be treated separately in 
Section V. 

When x' and x" are spatially separated, we may close the contour of integration over the upper half complex p-plain. 
Changing the integration variable to m = ip, Eq. (2.27) may be rewritten in the form. 



a{t')a{t")G+{x',x") = 
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where the contour of integration C on the complex u-plane is shown in Fig. 3. If we note the relations, 

P-,-{z')P:),{z") = P:),{z')P-,-{z"), P-,-{z')QUz") = P:.{z')Q-r{z"), (2.29) 

for integer n, we easily see that the residues of poles at u = n (n = 1, 2, 3, • • •) cancel out completely. Thus the poles 
which contribute to the integral are classified into three classes. 

class{a); poles of l/sin7r(u — v'): 

^-^ ^"^^ 1 n = 0,l,2,---. (0 < I/' < 1) ^^-^^^ 



class{b); poles of Q^,"{z"): 



class{c); poles of 1/ sin7r(u + u'): 



u = v' + n, n= 1,2,3,- ••. (2.31) 



The contribution from a class (a) pole is given by 

1 (-1)" 



an ■■= 



Stt sinh ( sin irv 
By using the following formulae 



-{u' + n)e-(^'+")f (^p-y+''\z')P;::+''{z") - ^ simru'e-''"''p-y+''\z')Ql'+''{z")^ . (2.33) 



Qu'^'^iz) = ^4 ^P:>{z); n = 1, 2, 3, • • • , 

^^""^^^ = l£^ (^^'""(") - ri = 0,1,2,. ■■, (2.34) 



we find 



87rsmn(; smTTi^' 

a„=0; n = 1,2,3,- ••. (2.35) 
The contribution from a class (b) pole is given by 



Snsmhc^ smTrz^' n sm27ri/' 

The limit in the square bracket is evaluated as 



lim {u-u' -n)Q-,''{z") 



Km (u-u' -n)Qz:'{z")= lim {u - {v' + n))V{v' - u+ ^ ' 



The contribution from a class (c) polo is given by 

-{v' - n)e("'-")^ {P-," +"{z')P;;, -"(^") - - smmy'e-''"''p-/ +"(^')<9;^'""(^")J • 



Stt sinh sin tti/' 



(2.36) 



M^i/'+n M^i/'+n r(j/' — U + 1) 

P;,'+"(2"). (2.37) 



(n- 1)! 2sin7r!/T(2z/' + n+l) 
Then we obtain 

<" ' - B.L( (n - l)!swtrnCrpl- + n + .f ^:*"^^')^'*"^''^' n = (2.38) 
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Using the second formula in Eq. ( 2.34| ), it reduces to 

^» = « - n)e(-'-")'^ ^^^'^^~" + ^^ p--'+"(z0^.r'+"(z"); n = 0, 1, 2, . . . . (2.40) 

Stt smh ^ sm TTi/' n! 

Combining these results, we obtain the Wightman function expressed in terms of the series sum of the residues. 
Here we focus on the case when — l/2<i/'<0or diiy' = —1/2. Then we have 

f oo oo 

^ ' 87ra(t')a(i")sinhCsin7ri/' 1 sin27ri^'^^^ ' ' J 

=:G+,)(x',x") + G+j)(x',x"), (2.41) 

where 



g -(i.'+n)C 



" (n - l)!r(2i/' + n + 1) ^ ^ ^ ^' 

S„ := ^'^^"'^'^ '"^17^"' ~ " ^ ^^ P,r'+"(zO^.';'^'+"(^'0, (2.42) 

and (a;', a;") and G^^(a;',x") correspond to the summations of yl„ and i?„ parts, respectively. 

In the next section, we shall show that the above expression indeed coincides with the Euclidean vacuum Wightman 
function. The case v' > Q will be discussed in Section IV. 



III. EUCLIDEAN VACUUM WIGHTMAN FUNCTION 



Equation (2.41) has been deriv ed by assuming that the Klein-Gordon inner product can be evaluated in terms of 
the expression given by Eq. (2.17). If this assumption is valid, we expect Eq. (2.41) to coincide with the well-known 
Wightman function for the Euclidean vacuum pO|] , 



G+(z',0;z",C) 



El 

87r2 



where 



u := s's" coshC 



1 = sinhi. 



(3.1) 



(3.2) 



Now we prove that Eq. (2.41) is identical to the Euclidean vacuum Wightman function. For convenience, we denote 



the function given by the series expression (2.41) by (x' , x"). The outline of the proof is as follows. First we show 
that {x' ,x") is de Sitter- invariant, i.e., it is a function of the invariant quantity u alone. Then, since G'^{x' ^x") 
is guaranteed to be a solution of the field equation with respect to both of the arguments x' and x" by construction, 
the form of G+(z', z", C,) is restricted to be a linear combination of two independent solutions as 



G+(z',z",C) = 



8^ 



aP;.\u)+P^—Q-}{u) 



(3.3) 



where a and /3 are some constants. Examining the asymptotic behaviour of the Wightman function for the limit of 
large spatial separation, a and [3 are determined as {a = l,/3 — 0}, which is of the Euclidean vacuum Wightman 
function. 

To show the de Sitter invariance of G+(z', z" , we use the fact that if a function F{z' , z", satisfies the relations. 



au9F(z',z",C) dudF{z',z"X) 



dC dz' 
du dF{z',z",C,) 



dz' dC 

du dF{z', z", C) 



dz" 



dz" 



dc 



(3.4) 
(3.5) 



F{z' , z" ,(^) is a function of u only, i.e., F(z' , z" , () is a de Sitter-invariant function. Of course we can directly check 
the relations (3.4) and ([3.E|) for F{z', z", C.) = G^{z' , z", C), but it is easier to show these relations first for 



8 



sin TTij' \ sin y-TTjy' 



\ n— 1 n— / 

Then if F{z' , z", C) is shown to be a function of u alone, G'^{z' , z" , C) can be obtained as 



s's" sinhC 

Hence G+(z', z", Q is also a funct ion of onl y u if z", C) is so. 
Now let us prove the relations (3.4) and (3.5). We define 



if(i) 



9m 

dudA^ s" 
du dB„ 



- cosh C - z" {y' + n)P;, +"(z') 



-(i/'+n)C 



(3.6) 



(3.7) 



^ sinhC ((i^' + n)z'P;;+"(z') + s'P;;+"+i(z')) 



(n- l)!r(2z/' + n+ 1) 

g-(i/+n)C 



(n~ l)!r(2z^' + ?i+l) 



p;;+"(z"), 



9z' 9C 
du dBn 
dz' s 



-f coshC - z" ) [u' - n)P,-;'^'+"(zO ^ ^^"''^'^ "^^^1^(2^^ " + ^^ P^r'+"(z'0: 



^ SinhC (-(^' - n)z'p-^'+"(z') + s'P;/+-+\z')) 



(_l)V(-/-«Kr(2i.--n + l) ,^ 



(3.8) 



Then 



2E(^i^^-^i.'') =E -2if e'; + 2z" (.' + ,,)p;'+"(,')-5"(e';-e-C)P;>+i(.') 

n=l n=l [ \ * / 



E 



(n- l)!r(2z^' + n + 1) 
1 



n!r(2j/' + n + 2) 



+2 



(n- l)!r(2t/' + ?i + l) 
1 



(^{v' +n+ i)fJp;;+"+i(z') + p;;+"+2(z')j s"p;;+"+i(z") 
p;,'+"(2')2"^;/+"(^") 



-P;,+"(2')s"-F'."+""'(^") 



-(i/'+«)C 



(n- 2)!r(2t/' + n) 

Each n term in the above expression is shown to be zero by iteratively using the recursion formula, 



v' 

In the same way, we have 



K +"+'(2) + 2(/^' + n + 1)-P;, +"+i(z) + n{2v' + n + 1)P;, +"(z) = 0. 



(3.9) 



(3.10) 



2E (4^^ - = 5: 2^^^^ - 2z" ) (.' - n)p-'^'+"(z') - ."(eC - e-';)P;;^'+"+i(z') 

n=0 n=0 \ * / 



(_l)ne(-''")Cr(2^'-n+l) „_,,^„^ 



E 



Y{2v' - n) 
(n + 1)! 



2(1/' - n - i)::^p-'^'+"+i(;2') - p-'^'+"+2(z') s"p;;+"+i(z") 



2 ('^' - n)r(2t/' - n + 1) p-y+n^^jy,p-y+n,j, 

Y(2v' - n + 2) D-i^'+n/ //p-i/'+ri-l/ // 



(n-1)! 



(-l)"e 



"p('^'-")C 



(3.11) 
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and each n term is shown to be zero by using the formula, 



P. 



(z) + 2{-v' + n + l)^P,;'^'+"+i(z) + [n + l){-2v' + n)P;y+''{z) = 0. 



(3.12) 



Thus we find that F{z' , 2", Q satisfies the condition (3^). The condition (3^) is also shown to be satisfied by replacing 
the roles of z' and z" . Hence we find that G'^{z' ^ z'\ Q is de Sitter-invariant as well as F{z\ z" , Q. 

As mentioned before, since G^{z' , z " , C) is de Sitter- invariant, the fact that this function satisfies the field equation 
restricts its form to that given by Eq. ( |3.3[ ). Moreover, as is apparent from the above proof of the de Sitter invariance, 
"^"i s-^^ both de Sitter-invariant and satisfy the field equation independently. This means 

that these two functions give two different de Sitter-invariant functions. To determine a and (3 for each of them, we 
examine their asymptotic behaviours at large C limit. In this limit, 2u s' s" e'' . 

The asymptotic behaviours of G^j(z', z", Q and G^^(2;', z", C) are determined by their leading terms. Noting that 



P. 



(^) 



-l/'-l 



r(i + z/') 



and 



r(2.) = ^rHr(z.-K- 



the asymptotic forms are explicitly written as 



27r3/2(iy' + l) 



47r(2u) -sm-Kv' sm2-Kv'V{v' + \)V{y' + 3/2)[r(-i^')]^ 



(l+"'')C(g'5")-(!^' + l) 



47r5/2 

471(2^) 
- J^L 



r(-!/'-l/2)r(i.' + 2)(2w)-'^'-2, 

.T(.- + l/2) vc. 
0Fsin7ri^T(j/' + 1) ^ ^ 

r(i.' + l/2)r(-!y' + l)(2uf'-i. 



Comparing these with the asymptotic form of the right hand side of Eq. (3.3), 

r(j/' + 2)r(-j/' + i) ^ , /3 



aP-Mu) 



47r5/2 



ar{i'' + l/2)T{-iy' + l)(2u)'^'-^ + (a + /3)r(-i^' - l/2)r(iy' + 2){2u) 



-u'-2 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



it is easy to see that G^^(z',z",C) corresponds to the choice of the coefficients {a = 1, (3 ~ —1}, and G^j(z',z",C) 
to {a = 0,/3 = 1}. Thus the sum G+ — G^^ -I- G^^ correctly gives the Euclidean vacuum Wightman function. This 
completes our proof. 



IV. EXTENSION TO > 



Although we have assumed —1/2 < v' <Q or '^v' = —1/2 to derive the expression (2.41) for the Wightman function, 
if we examine the proof in the previous section that it coincides with the Euclidean vacuum Wightman function, we 
find it is valid irrespective of the value of v' , that is G+ = G^ as long as the series converges. Then if we consider 
the case v' > 0, we find that G+(z',z",C) differs from the o riginal G"''(z', z", C), whi ch is given by the integral of 
products of the mode functions over real values of p, Eq. ( 2.27 ), or the contour integral ( 2.28 ), because different poles 
contribute to the integral. Specifically the contribution of ao comes in and cq goes out as v' becomes positive. Thus 



G+(z', z", C) - G+(z', z", C) = (ao - co) 

s's" 
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Sns's" sinh ( 



tll,'<p-/(,')p-/{,") 



-4.5/.r(-' + W + i/2)— ^(.vr-. (4.1) 

Phys ical o rigin of this difference between v' < and > may be understood as follows. The inverse transformation 
of Eq. ( |2.25| ) is given by 

±1 P Jq smhC 



(T=± 



As seen from the asymptotic behaviom' of the Euclidean vacuum Wightman function for large C given in Eq. (3.17) 
with {a — 1,(3 = 0}, the integral in Eq. (4.2) becomes manifestly divergent for i/' > if we set G"*" = G^. This 



means that the Euclidean vacuum state cannot be described by any power spectrum in the usual sense in which the 
summation of products of the mode functions is taken over real values of p. In deriving the normalized set of mode 
functions in Section III, we did not assumed anything other than that the Klein-Gordon norm can be evaluated on 
two isolated hypersurfaces denoted by (I) and (II). Thus we conclude that the discrepancy at > is due to the 
breakdown of this assumption. 

Let us consider a linear functional space, J-', which consists of all normalizablc functions which satisfy the field 
equation. Here the normalizable functions mean that they have finite Klein-Gordon norms, apart from the trivial 
factor of the (5-function. The Klein-Gordon inner product is defined on a Cauchy surface and is independent of the 
choice of the surface. To evaluate it, there is one convenient choice of the hypersurface, (IV), which is specified by 
rc = constant and included in the region C as shown in Fig. 2. 

As the system has the 0(3, 1) symmetry, it is natural to take the basis vectors of !F from functions of the separated 



form as given in Eq. (2.13). Restricting to be positive, if we evaluate these mode functions on the hypersurface 
(IV), they diverge at tc ~ 7r/2 or tc = —tt/2 with infinite cycles of oscillation. However the oscillatory divergence is 
not so seriously bad that the Klein-Gordon norm still has a finite value. It then follows that two modes which have 
different separation constants p^ and p'^ are orthogonal. Therefore they are contained in J- and may be chosen to be 
the basis vectors. Moreover, as the mode function fpi(rc) decreases rapidly as tq increases, after the deformation of 
the hypersurface (IV) to that consists of (I), (II) and (III), the contribution to the Klein-Gordon norm fro m (II I) can 



be neglected as was assumed before. Thus the Klein-Gordon norm is correctly given by the expression ( 2.16 ). For 
completeness, in Appendix A we show that the Klein-Gordon inner products evaluated on the hypersurface (IV) are 
finite for mode functions with p^ > Q and in fact coincide with those evaluated on (I) and (II) . 

On the other hand, for p^ < 0, almost all modes diverge at ^c" = tt/2 oi tc = —tt/2 faster than {t ^ tt/2)^^ 
(the critical power exponent depends on the dimension) and their Klein-Gordon norms diverge. Hence they are not 
contained in J^. However, when v' > 0, there exists an exceptionally set of modes with p = iv' which behave well 
enough, 

, P_\y,']jisivL\\rc) 

U(*)im = H{-icostcY -1 Zi^^ Yi„M. (4.3) 

V « cosh rc 

Its Klein- Gordon norm is finite and is given by 



icosh^rc f"^"^ ,^ ^ f ir^ ( 9^i*)im duu)i,,n' 

' dtc cos^- I '^<> ^ — — — -- 

/2 

20Fr(j/'; 



(u(*)/m, -Jp J dtcCOStc y '^^l g^^ U(*)l'm' ~ U{*)lm 



" T{v' + l/2)T{-v' + / + l)r(z/' + / + 1) 

= : N(^^)i5iv5^,n' ■ (4.4) 
In the region i?, these mode functions take the form, 

/ , -P_,t_V/9(coshr/j) 

u^^)i„Xx) = Hisinhtj^r )!L : YUn). (4.5) 

Vsmh rji 
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They have exactly the same functional form in the region L as well, with {tn^rn) replaced by {tL,rL)- Therefore 
the contribution to the Wightman function from these modes in the region R or L becomes (with the suffix R ot L 
suppressed) 



G+)(z',z",C) :=E 



N, 



l,m 

, p-'rY',(coshr') p-'rV/o (cosh r") 



8^3/2 



r(-zy' + i)r(i/' + 1/2) zy' 



V sinh r' \J sinh r" 

,Fr;/_\/,(coshc) 



-\/sinh C 



^-J/-i/2(coshr) 



lim 



\J sinh ? 



(s's")'''-' 



r(-i.' + i)r(i/' + 1/2) 5^^^ (s'5")^'-\ 



47r5/2 * ' - ' ' ^' ' sinhC 

where, in the third equality, the relation required by the 0(3, 1) invariance, 

Y{-v' + / + \)Y[v' + Z + 1) -Pr^rii^/Vcoshr') P~_'jX^(cos\^r") 



(4.6) 



l,7n 



r(-j/' + i)r(i^' + 1) 



V sinh r' 

1/2 



1/2 

V sinh r" 



Yimin')Yim{n") 



P:'J' [cosh 



\/sinh C 



lim 



-P- ^^''^-1/2 (cosh r) 
Vsinh r 



2"^r.V-i/2(coshc) 



TT Vsmh^ 



(4.7) 



has been used. For i/' — —ip {p > 0), this is just the equation describing the completeness of Ypim, Eq- ( 2.26 ), which 
is well-known to hold. In Appendix B, we show it holds for i^' > 0, ly' ^ integer as well. Thus we find that 



G+(z', z", C) - G+(z', z", C) + G+) (z', z", C). 



(4.8) 



As clear from the above, the degrees of freedom described by cannot be quantized in the open chart, where 

these modes are not normalizable. However if the quantization is performed on some other hypersurfaces, say, the 
time constant hypersurfaces of a closed chart, we do not confront with any problem in the quantization. We may then 
quantize the field there, but instead of using the harmonics on the three-sphere, we may expand the field operator in 
terms of the mode functions associated with an open chart. 



(x) = / dp y2 Vpcrlm{x)dpalm +y2'^{*)lma(*)lm + h.C; I > 1^' > Q, 

JO , I 

(7,l.ra Lm 



(4.9) 



where the orthonornialized mode functions are denoted by ua as before (see Eq. (2.24) for v' < 0) and are now given 
by 

Vpalrnix) = /p; (r) F;™ (fi) , 



y^N^ sinhi 



H , P \ ^//o(coshr) 

(sinht)-'-^ -"-y'] 'yU^). 



\J sinh •) 



(4.10) 



In the above, the coordinates (i, r) stand for {tii,rii) or {tL-,rL) if x is in the region R or L, respectively, and for 
{i{±tc — 7''/2), rc + 17^/2)) in the regi on G where the upper (lower) sign corresponds to the analytic continuation from 
the region R (L) as specified in Eq. (2.3). 

Thus we have found a complete prescription to express the Euclidean vacuum state in terms of the mode functions 
associated with an open chart. Even in general situations in which the background differs from the exact de Sitter 
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space and the mass of the scalar field changes in time, the above method of decomposing the field operator in terms of 
the mode functions which respect the 0(3, 1) symmetry will be a very useful tool to evaluate the evolution of the field 
as long as the quantum state is initially in the Euclidean vacuum. Furthermore, we expect that the logical procedure 
given here should play an important role when we consider various problems associated with the quantum tunneling 
through 0(4)-symmetric bubble nucleation. 



V. CONFORMAL AND MINIMAL MASSLESS CASE 



In this section we consider the cases i/' = and i'' ^ 1, which correspond to the conformal massless case and the 
minimal massless limit, respectively. 



A. Massless conformal coupling case 



First we consider the case ly' = 0. In this case the integration of Eqs. (2.25) and (4.2) can be easily performed 
because the Legendre function reduces to an elementary function as 

where ij is related to t by the relation, sinhr; = — 1/sinht. Then, the right hand side of Eq. ( ^.25| ) may be expressed 



as 



7r2 J„ ^sinhC 1 1 - e'-^^P 1 - e-^^P 



^ / dp—B { ^—^e-'"'' " ' + . ^-.^y ' (5-2) 

This integral can be easily done by closing the contour and evaluating its residues. Then we find that it coincides 
with the Euclidean vacuum Wightman function, 

G+(z',z",C) = ^T^. (5.3) 

If we focus on either of the region i? or L, a natural choice of the positive frequency function is Xp = e.~^P^ / ^/2p. 
Interpreting Eq. ( |5.2| ) in terms of these mode functions, the Euclidean vacuum is seen as a thermal state, as was 
pointed out in pl| . This thermal nature may be understood as a consequence of the loss of information. For complete 
description of a quantum state, it is necessary to specify the states both in the regions R and L. Hence if we 
consider the expectation value of an observable which has support only in the region R (or L), it does not depend 
on the quantum state of the region L (or R). This effectively brings a pure state into a mixed state. Of course, 
this phenomenon is not peculiar only for i^' ~ but it appears in the easiest form for this special case. In fact, the 



expression for Xpa given in Eq.(2.21) does already indicate the thermal nature; the ratio of the squared absolute values 
of the coefficients in front of ("positive frequency" in the natural sense) and P^^,^^ ("negative frequency"), after 
normalized by Npcr and summed over cr = ±, is e^'^P. 

B. Massless minimal coupling limit 

Next let us consider the massless minimal coupling limit, v' 1. Using the following formulae, 

r(3 + e) = 2 + e(3-27) + 0(e2), 

r(-e) = -i-7 + 0(e), (5.4) 

for e — > 0, the v' 1 hmit of G^{z' , z" , C) is given by 
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1 



87r2 1 z/' - 1 



1 



u+1 



-In 



u+1 



- 2 



+ 0(^'-l) . 



In the same limit, the contribution from p — iv' modes to the Wightman function given by Eq. ( |4.6D becomes 

1 



87r2 \v' -I 



C coth C + In 



G+)(z',z",C)- 

On the other hand, for v' = 1, Eq. ( 2.25| ) reduces to a very simple form 

1 



+ 0{u' -I) 



sinpC e^P jz' + ip)(z" - ip) 



sinh C sinh vrp I + 

The above integration can be done elementarily to give 

'1 + It 



G+(z',z",C) = ^{T^-ln 
Svr^ L 1 + M 



In 



+Ccothc}. 



(5.5) 



(5.6) 



(5.7) 



(5.8) 



Adding Eqs. (5.6) and (|5.8|), we immediately see that G+(z',z",C) — G%{z' ,z" X)- We note that the Wightman 



function same as Eq. ( |5.8| ) has been obtained by Allen and Caldwell ^] in a totally different fashion 



VI. CONCLUSION 



We have derived the expression for the Euclidean vacuum Wightman function of a scalar field in de Sitter space 
in terms of the mode functions associated with the harmonic functions on an open chart, Ypi^m i.e., on hyperbolic 
time slices which respect the 0(3, 1) symmetry. The formula we have obtained is applicable to any mass or curvature 
coupling, provided the effective mass square is positive, i.e., M'^j^ = + 12^H^ > 0. Usually the power spectrum 
of the quantum fluctuations is described by the modes with real values of comoving wavenumber p. However we have 
found that the Euclidean vacuum Wightman function cannot be described in this standard manner for Af^^j < 2H^. 
We have found that an additional contribution from the modes with an imaginary value of p is necessary to describe 
the Euclidean vacuum Wightman function in this case. These modes are not square-integrable on a hyperbolic time 
slice in an open chart but their Klein-Gordon norms are finite if we take a complete Cauchy surface in de Sitter space 
to evaluate them. 

The expressions for the Euclidean vacuum Wightman function and the field operator we have obtained will be a 
powerful tool to calculate observable quantities in open universe inflation models such as the cosmic temperature 
fluctuations when the initial quantum state is in the Euclidean vacuum [p|-|To|. Furthermore, the technique developed 
in this paper to obtain a set of orthonormalized mode functions associated with an open chart should be applicable 
to more general situations provided the system has the 0(3, 1) symmetry. The issue of the quantum state after false 
vacuum decay through 0(4)-symmetric bubble nucleation is one interesting example of such situations. This issue 
has been discussed and investigated by several authors 1 13 17| but in a rather formal manner. In particular, although 
a method to determine the unnormalized mode functions after tunneling has been given, the normalization procedure 
applicable to general situations has not been explicitly given, except for the case of flat spacetime background . 
Application of our normalization technique to some specific models of false vacuum decay will be discussed in a 
forthcoming paper. 
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APPENDIX A: 



Here we show the Klein-Gordon inner products on the hypersurface (IV) are finite for the mode functions with 
> and they coincides with those evaluated on the hypersurfaces (I) and (II). 
We consider the mode function Upo-Zm given by 



a(tc) 



(Al) 



where a(ic) — H ^cosic, Xp,o- is giv en by Eg . ( 2.19| ) and f pi is by Eq. (2.11), both with relevant analytic continuations 
to the region C specified in Eq. (2.3). In particular, the analytic continuation of fpi to the region C is given by 



fpi{rc) 



r(zp + / + 1) P^p-l|2^^ sinh rc) 



\/ i cosh rc 



(A2) 



which plays the role of the positive frequency function. The Klein-Gordon inner products among Upaim on the 
rc = constant hypersurface (IV) take the form. 



palm ^p' a' I'm' 



'^palm J ^p'cr'l'm' 



Tr/2 



i cosh rc 
IP 

Kl^Kllp'<r'Su'Smm' , 



dtc cos tc / dil 



-7r/2 



palm 



Ore 



Up' (j' I'm' '^palm 



d Up' a'V m' 

Ore 



(A3) 



where 



N. 



.(1) . ,2 (9fpi-j^ f dfpi\ 
pi :=*cosh rc^-^fpi-fpij^\ 



drc J 



N. 



(2) 



p(T,p'(T' 



COS tc 



Xp,aXp' ,a' 



(A4) 



(2) 

In the above, we have anticipated that the factor N^J should contain the delta function 5{p — p') so that we may 



put p' — p when evaluating N^^ . This will be justified shortly. 



The evaluation of the factor N^'^' is straightforward. We immediately obtain 



(1) _ 2p 



pi 



sinh np . 



(A5) 



The evaluation of the factor N^J ^i^, needs a bit of consideration. First, from the symmetry that Xp,a{ic) = Xp,a{<^tc) 

f 21 

(ct = ±), we see that N^J^,^, must be proportional to 5aa'- Next, by using the field equation for Xp,(t, 



— cos rc^T 

otc otc cos rc 



1 / 9 9 , 9 

+ {v^ - -)cos^ rc 



Xp,<y = (1 +P^)Xp,'T 



(A6) 



we find 



Tr/2 



(1+p^)/ dtcicOStc) ^Xp.'yXp',a' = + P'^) I dtcicOStc) ^Xp,crXp\ 



Tr/2 

r(2) 



tt/2 



Tr/2 



which implies that N^J^,^, = if p' ^ p^- Hence it should be expressed as 

r(2) 



(A7) 



(A8) 



where we have restricted p and p' to be positive (negative p modes are equivalent to positive ones because of the 

f 2) 

symmetry X-p,cr = Xp.a)- Thus our remaining task is to evaluate Nj 



(2) 

It turns out that Npa for real p can be evaluated without detailed knowledge of the behavior of Xp,a because the 
divergent contribution a.t p' = p which determines the coefficient of the J-function comes only from the boundaries of 



per 
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integration at tc = ±7r/2. Therefore, to obtai n Np^ , we only have to know the behaviour of Xp,cr near the boundaries. 
For this purpose, we use the expression (2.21) for Xp,aj which we recapitulate: 



2 sinh irp \ T^v' + ip + 1) 



Xp,cr — < 



— P'^(zr) 

2 sinh irp \.r{iy' + ip+ I) " ^ ' r{iy'-ip+l) 
Then noting the asymptotic behavior of near the boundaries, 

we find 



r{v' -ip+1) 



r(i - ip) 



r(i - ip) 



r(i - ip) 



r(i - ip) 



Xp,<y 



api-tc + n/2)-''P-l3pi-tc + 7r/2)'P for tc ^Tr/2-0, 
a{ap{tc+7r/2)-'P-Pp{tc + TT/2yP) for tc ^~7r/2 + 0, 



where 



-■n-p/2( TTp 



{e^P - ae-''^'' ) 



/9p 



2sinh7rpr(i^' + ?p+ l)r(l - ip) ' 
2sinh7rpr(z/' - + l)r(l + ip) ' 



Then 



Nj,l^Sip-p') = lim 



dx 



X 



dx 

X 

|2 I \fD\2 







from which we obtain 



2n{\ap\^ + \Pp\^}6ip~p') 



(2) 2 (cosh np — a cos v'lr) 
P" " psinh7rp|r(i/' + 1)|2" 



From Eqs. (A5) and (A14), we finally find 



4(cosh7rp— (tcos^^'tt) , 

J^pcrlm,p'a'l'm' = i-r^^ , . ■ r~rTl9 "[P-P )0aa'Owd„ 

n \l [I'' + ip + Ijp 



We see that this exactly coincides with the result obtained in Eq. (2.22). 



(A9) 



(AlO) 



(All) 



(A12) 



(A13) 



(A14) 



(A15) 



APPENDIX B: 

In this Appendix we show the outline of a proof of Eq. ( [4.7D . After summing over m in the left hand side of the 
equation, the relation to be proved becomes 

- Pr;/_\ (cosh C) P:'J\^,icoshr) f2P:'J\/,{coshC) 

y Ml = -p= hm = \\ -p= , 

^ VsinhC '■^0 Vsinhr V tt VsinhC 

M ri-.' + l + m.' + l + l)P:l7y,%{coshr')P:l7X%{coshrn 

■ r(-.' + l)r(.' + 1) Vsinhr" (2^ + l)^'(cosO), (Bl) 
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where 

cos e = cos e' cos e" + sm e' sin e" cos(0' - 0"), 

cosh C = cosh r' cosh r" — smh r' sinh r" cos 6, (B2) 

and Pi(cos8) is the Legendre function. Here we assume v' > Q and v' ^ integer. 

The strategy is similar to the technique employed in Section III. First we show that Mi converges. Then we 
show the result depends only on C. Then taking the limit r" 0, it is readily seen that Eq. (Bl) holds. 

First let us show the convergence of J^i ^'^i- For this purpose, we rewrite Mi as 

where the function Fi{r) is defined by 

V{-u' + 1 + 1) PZl7l[^/^{coshr) 



Fi{r) 



T{-v' + l) Vih^ 
2 



T{-v' + l)sinhye''^'' Q-" (coshy), (B4) 

TT 



and y is related to r as sinhy = 1/sinhr or y = lncothr/2. With the aid of the integral representation, 

- / / ^/TTsinh"''' y f°° e^^'+^/^^^da; 



V2T{v' + 1/2) Jy (cosh a: - cosh ?/) -'"'+1/2 ' 
for W > -1/2, 3?(i^'-0<l, 2/>0, (B5) 

we see that F;(r) is real and positive for Z > ly' — 1. Further it forms a monotonically decreasing sequence for Z > v' ~1. 
Hence it is bounded as 

< Fi{r) < e-('-")^f^„(r) = (tanh i^„(r), (B6) 
for I > n — [v'] where [v'] is the largest integer not exceeding v' . Also for the same n, 

ri-Xtll) < + + ^ - 1) • • • + ^ - 2n - 1) < + lf-+\ (B7) 
Thus, for sufficiently large I, the absolute value of Mi is bounded from above as 

\Mi\ < + (^tanh^tanh^y""K(r')K(r"), (B8) 

which becomes exponentially small for large I. Therefore the series J^i manifestly convergent. In particular, if 

Q < v' < \, Fn{r) is bounded from above for real positive r and approaches zero as r ^ 00. Thus the series converges 
uniformly and absolutely for < i^' < 1. Furthermore, even for > 1, as long as we keep r finite, say r < r^axi 
the series also converges uniformly and absolutely there. In this sense, the convergence is uniform for any v' > 0, 
1^' integer (uniform convergence on compact sets). 

Next we show the series J2i indeed a function of C alone. The equation to be proved is 

V-/ 2^ ^f dC, dMi dC dMi 1 „ 

^ cos^9-l \ T^-TT^T^ in\=^- ^9 

^ a cosh r' a cos B a cos B o cosh r' J 

This calculation is a bit more complicated than the previous case in Section III because the counterpart of e'^'^'^^"^'' 
there is Pi(cos0) here. Taking derivatives with respect to cos 8 gives terms as cos BP; (cos 9) and cos^ BP; (cos B). 
The factor cos B must be removed from these terms by using the relation, 

cosBPi(cosB) = ^^{(/ + l)P;+i(cosB) + ZP,_i(cosB)}. (BIO) 
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After carrying out this procedure, the left hand side of Eq. (B9) takes the form, J^i Ci{r' , r")P/ (cos Q). Then applying 
the recursion relation for the Associated Legendre functions, it can be shown that Ci{r' , r") = 0. Thus the right hand 
side of Eq. (Bl) has the desired invariance, which then immediately implies it is equivalent to the left hand side. 
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FIGURE CAPTIONS 



Fig. 1. Conformal diagram of de Sitter space with the coordinates (i,r) = (ti?,r-fl), {IltTi,) and {tc,rc) spanning the 
regions R, L and C, respectively. The thin real lines denote the surfaces {r = constant} and the broken lines 
{t = constant}. 

Fig. 2. Two different Cauchy surfaces on which the Klein-Gordon innner products are evaluated. One of them consists 
of hyperbolic time slices (I) and (II) in the regions R and L, respectively, and the hypersurface (III) which 
connects them. The other one (IV) is an rc = constant hypersurface entirely contained in the region C. 

Fig. 3. The contour of integration on the complex u-plane for Eq. ( ^.28| ). 
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